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Abstract 

The problem of persistence of four-frequency tori is considered in 
models represented by the coupled periodically driven self-oscillators. 
We show that the adding the third oscillator gives rise to destruction of 
the three-frequency tori, with appearance of regions of either chaotic 
attractors or four-frequency tori. As the coupling strength decreases, 
the four-frequency tori dominate, and the amplitude threshold of their 
occurrence vanishes. Also, for three oscillators, a domain of complete 
synchronization of the system by the external driving can disappear. 

PACS: 05.45.-a, 05.45.Xt 



1 INTRODUCTION 

A problem of effect of external driving in systems of mutually coupled self- 
oscillators is a fundamental part of the theory of synchronization important 
for a variety of applications [1], |2]. It is a part of a field related to birth and 
synchronization of multi-frequency quasi-periodic oscillations. In particular, 
it goes back to the concept of the onset of turbulence advanced by Landau 
and Hopf, which is based on successive appearance of growing number of 
quasi-periodic components of the motion with increase of the Reynolds num- 
ber [3]. As known, further developments of nonlinear dynamics and bifurca- 
tion theory have made this hypothesis doubtful in some respects. According 
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to Ruelle and Takens [1], in systems with multiple oscillatory components 
arising successively due to a sequence of bifurcations, the quasi-periodic mo- 
tions associated with torus attractors cannot occur as typical at a number of 
frequency components larger than three; instead a strange attractor should 
appear. On the other hand, computations for model systems undertaken by 
several authors do not quite agree with this statement: in certain situations 
tori of sufficiently high dimension seem persisting and observable [5]. Thus, 
the question of how the transition to chaos in multidimensional systems de- 
velops, and what is the role of the number of oscillatory components, or 
degrees of freedom, involved in the dynamics, deserve further studies. The 
phase approximation whereby dynamics of coupled periodically driven self- 
oscillators is investigated appears to be effective and illuminating in the study 
of this problem [5]. In this paper, in the framework of the phase equations 
we discuss transformation of the regions of the different regimes in the pa- 
rameter space when the number of driven oscillators increases. Thus, we will 
be interested in the question of how the transition to high- dimensional tori 
and chaos depends on the number of oscillators in the system. Also, we show 
that an increase in the number of oscillators up to three leads to a new effect 
- a domain of complete synchronization may disappear. 

2 PHASE EQUATIONS 

Let us consider a system of three dissipative coupled van der Pol oscillators 
driven by an external periodic force (Fig.l): 

X — (A — x'^)x + X + fi{x — y) = B sin ut, 

y-{X- y^)y + (1 + Ai)y + fi{y - x) + fi{y - i) = 0, (1) 

z-{X- z^)z + (1 + A2)2: + ii{z -y) = 0. 

Here A is a parameter responsible for the Andronov - Hopf bifurcation in 
the independent oscillators, Ai 2 are frequency detuning of the second and 
third oscillators from the first oscillator, fi is the coupling coefficient, B is an 
amplitude of the external periodic force, u is its frequency. The frequency of 
the first oscillator is assumed to be normalized to unity. 

The system ([1]) may be analyzed in terms of complex amplitudes (quasi- 
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Figure 1: Schematic representation of a system of three dissipative coupled 
oscillators driven by an external periodic force. 



harmonic approximation) if the parameter A is small. For this, we set 

z = c(t)e*^* + c*(t)e-^'^*. 



(2) 



Here, a, b, c are the complex amplitudes of the oscillators slow varying in 
time scale of the basic oscillations with unit frequency. We use the additional 
conditions traditional for this method 



0, 



(3) 



Then, we substitute ([2]) and ([3]) in ([T]), multiply the equations by factor 
and perform averaging over a period of the basic oscillations. Assuming 



-iujt 



w ~ 1 the resulting equations for complex amplitudes are as follows 

b 



4' 



2a = Xa — \a\ a — 2i^la — fi{a — b) — 

2b= Xb- \b\% - i{2n - Ai)b - fi{b - a) - fiib - c) 
2c = Ac - |cpc - i{2n - A2)c - ^(c - b), 



(4) 



where = a; — 1 is detuning parameter of the external driving force with 
respect to the first oscillator. Note that in the equations dl]) the parameter 
A may be eliminated via variable change, and then the other parameters will 
be renormalized with the factor A~-^ [5]. Let us set a = Re^^^, b = re^^^, 
c = ve^^^, where R, r, v are real amplitudes, and "012 3 are phases of the 
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oscillators. Following the Refs. [T] , [1] , [S] , we assume that the motion takes 
place close to the unperturbed limit cycles of the oscillators, and set i? ~ 1, 
r ^ 1, V ^ 1. Then, the phases have to satisfy the equations 

■01 = - ^/isin('!/)i - 1P2) + bsinipi, 

11 1 
= -^^ + + 2^ sin(^i - ^2) - sin(^2 - ^3), (5) 

ip^= -Q + ^Aa + ^/isin(^2 - ^3)- 

Further we shall use the relative phases of the oscillators 6 = ipi—il)2 and 
V9 = ■02 — "03- Equations (E]) generalize the equations obtained in Ref.|l] to 
the case of three driven oscillators. 



3 SYNCHRONIZATION OF TWO- 
FREQUENCY OSCILLATIONS 

First, let us turn to the case of two oscillators. Then phase equation are 

■01 = —VL /i sin 9 + h sin ipi , 

1 1 
■02 = -fi + -Ai + -//sin 6*. 

We fix the "internal" parameters of the system Ai = 0.3 and = 0.4. 
Obviously, in absence of the external driving 6 = 0, the condition of mu- 
tual synchronization of the subsystems is /i > [1]. Hence, the above 

selected parameter values correspond to the regime of the mode locking for 
the unforced system. The external signal applied can destroy this regime. 

To visualize the domains of different dynamical regimes the results of 
parameter scan are represented in Fig. 2a for the h) parameter plane. 
Shown are the region of complete synchronization of the phases of oscillators 
by the external driving P and regions of the quasiperiodical regimes, namely, 
two-frequency tori T2 and three- frequency tori T3. For each parameter value 
we perform computations to reach a sustained dynamical regime of the model 
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([HD, and then evaluate two Lyapunov exponents, Ai and A2, using a standard 
Bennettin algorithm. The nature of the regimes is determined as follows. The 
case Ai < 0, A2 < corresponds to a complete phase locking of the oscillator 
phases by the external force. The cases Ai = 0, A2 < and Ai = A2 = C0 
are interpreted as associated with quasiperiodic regimes of two-frequency tori 
and three-frequency tori, respectively. 

It should be noted, that numerical tests are inevitably crude. In par- 
ticular, some of the points labelled as three-frequency quasiperiodic may 
well be high-order resonance two-frequency quasiperiodic, and four-frequency 
quasiperiodic may well be three-frequency quasiperiodic with high-order re- 
lation. Furthermore, such a Lyapunov classification labels as quasiperiodic 
some points, at which the Lyapunov exponent may vanish, for example, bi- 
furcation points, but this situations usually take place only at the borders of 
domains of stable dynamical regimes and do not influence the overall picture. 

The two-frequency regimes T2 correspond to a partial locking of the os- 
cillatory system by the external periodic force, while the relative phase os- 
cillates near some stationary value. Observe that these regimes dominate at 
small driving amplitudes, that is, they occupy a largest area on the param- 
eter plane there. At large amplitudes, again the two-frequency tori occur, 
but now the first oscillator is partially locked. In the intermediate domain 
of moderate amplitudes, the three-frequency tori T3 are observed, and the 
region of their occurrence is pierced by narrow tongues of higher order reso- 
nance two-frequency tori. 

4 SYNCHRONIZATION OF THREE- 
FREQUENCY OSCILLATIONS 

Now, let us add the third oscillator into the model system, and turn to the 
complete set of equations (|5]). We keep the same parameter values Ai = 0.3 
and = 0.4, and set additionally A2 = 1. In the absence of the external force, 
presence of the third oscillator gives rise to regime of beats at these parameter 
values. In other words, the third oscillator destroys the mutual mode-locking 
of the first and second oscillators. As a result, essential complication of the 
arrangement of domains in the parameter plane [Q, b) is observed. 

^Checking the zero Lyapunov exponents was carried out with an accuracy of 10"'^. The 
accuracy of calculating the main Lyapunov exponents was at least 10~^. 
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Fig. 2b shows the parameter plane chart for the system of three coupled 
driven oscillators. Now, the picture is drawn using numerical analysis of 
three Lyapunov exponents. New formations are observed on the diagram. 
These are domains of four- frequency tori T4, where Ai = A2 = A3 = (up to 
numerical errors), and domains of chaos C, with one positive exponent, that 
is Ai > 0, A2 < 0, A3 < 0. 

One can see that the domain of complete locking of all oscillators by 
the external driving force P decreases notably in width. Besides, a finite 
amplitude threshold of the complete synchronization appears. 

The regions of three-frequency tori disappear, and instead of them the 
four-frequency tori and chaos arise. Besides, a finite amplitude threshold of 
the complete synchronization appears. At small amplitude three-frequency 
tori replace the two-frequency quasiperiodicity. 

It is interesting to discuss the picture in the context of the Ruelle - Takens 
hypothesis. Adding the third oscillator obviously facilitates destruction of 
the three-frequency tori, but it can result in appearance either of chaos, 
or of the four-frequency tori. The domains of occurrence of these types of 
attractors coexist on the parameter plane. Both quasiperiodicity and chaos 
look typical, no one prevails. So, we do not observe an extrusion of the tori 
by the chaotic regimes, as one could expect, perhaps, basing on the Ruelle - 
Takens assertions. 

The winding-number chart of two-frequency tori shown in Fig. 3 demon- 
strates the complexity of the observed picture. 

In terms of phase equations ([5]), a two- frequency torus corresponds to 
a closed invariant curve in the phase box (-01, 1^2, ips), with boundaries (0, 
2n). The parameter scan in Fig. 3 is obtained as follows. For each parameter 
value we determine via numerical simulations a set of numbers p, q, r, which 
are essentially the numbers of intersections of the invariant curve with the 
edge surfaces of the phase box (with exclusion of returns). These numbers 
determine what we call the winding number of the invariant curve: w = p : 
q : r. Then, each point is attributed with a definite color depending on the 
winding number we got. The light gray tone designates regimes, which do 
not have a definite quantifier in the frame of the procedure, like the three- 
frequency tori, four-frequency tori, or chaos. The domains of main resonances 
are marked with subscriptions in the respective areas of the chart in Fig. 3. 
Note that the regimes of type w = : q : r correspond to partial locking 
of the first oscillator by the external force, and the regime w = 1 : 1 : 1 
corresponds to synchronization of all oscillators as the unified system. 
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Figure 2: Charts of Lyapunov exponents for the model systems of two (a) and 
three (b and c) driven oscillators. Parameter values are Ai = 0.3, A2 = 1, 
/X = 0.4 (a), fi = 0.4 (b), /U = 015 (c). The color coding is indicated in the 
bottom part of the figure. The marks fl2 and ^3 indicate frequencies of the 
second and the third oscillators respectively. 

Now let us consider transformation of the picture under decrease of the 
coupling between the oscillators up to the value /i = 0.15. The respective 
diagram is shown in Fig. 2c. In this case, the four-frequency tori dominate, 
and they effectively extrude the chaotic regimes, which occupy now only some 
extremely narrow strip-like areas placed mainly along the borders of three- 
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Figure 3: A chart of two-frequency tori for the system ([5]). Parameters are 
corresponded to Figs, lb and c. 

frequency tori domains. The four-frequency tori dominate as well at small 
amplitudes, so, the amplitude threshold for them vanishes. Additionally, 
there is a small chaotic region in a neighborhood of the frequency of the 
third oscillator, which has, however, a non-zero amplitude threshold. 

A characteristic feature of the Fig. 2c is disappearance of the regime of the 
complete synchronization P. Instead, a narrow strip of the two- frequency tori 
remains. The winding-number chart for this case shown in Fig. 3b illustrates 
this effect. One can see the regimes with winding number w = : : 1 
instead of the region of complete synchronization. This regime corresponds 
to partial locking of the first and second oscillator by the external force. 
This phase of the third oscillator demonstrates unlimited growth. Thus the 
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system of three oscillators can be configured so that complete synchronization 
is impossible. This is important difference from the case of two oscillators. 

In this paper we have considered the chain of three oscillators with exter- 
nal force effecting the outermost one. Otherwise, there are many possibilities 
for couphng and forcing: we may take the ring of oscillators or affect the oscil- 
lators simultaneously, etc. It may be conjectured that 4D-tori would persist 
for either configuration, but what concerns the possibility of absence of com- 
plete synchronization, the answer may depend on concrete configuration of 
coupling and external forcing. 

5 CONCLUSION 

A passage from two-mode to three-mode self-oscillatory system driven by 
periodic external signal gives rise to destruction of three-frequency tori with 
formation of regions of four-frequency tori and chaos. At moderate coupling 
parameter values neither of these regimes dominates. Decrease of coupling 
leads to almost complete extrusion of chaos by the four-frequency quasiperi- 
odic regimes. Besides, such regimes dominate at small amplitudes of external 
driving. New effects for the case of three oscillators is that at sufficiently small 
coupling the regime of the complete synchronization of three oscillators by 
the external driving disappear. 

The work has been supported by RFBR, grant No 09-02-00426. 
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